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Starting from few simple examples we have proposed a general method for finding an exact analytical 
solution for the two state scattering problem in presence of a delta function coupling. We have also 
extended our model to deal with general one dimensional multi-channel scattering problems. 



" ' INTRODUCTION: 

o : 

CS| , Nonadiabatic transition due to potential curve crossing is one of the most important mechanisms to effectively induce 

^ electronic transitions in collisions [H . This is a very interdisciplinary concept and appears in various fields of physics 

D . and chemistry and even in biology 0-0] ■ The theory of non-adiabatic transitions dates back to 1932, when the 

pLn ' pioneering works for curve-crossing and non-crossing were published by Landau 0, Zener ^ and Stueckelberg [13] 

, and by Rosen and Zener [ll[ respectively. Osherov and Voronin solved the case where two diabatic potentials are 



o 



fS| constant with exponential coupling 12|. C. Zhu solved the case where two diabatic potentials are exponential with 
exponential coupling [3]. In this paper we consider the case of two or more diabatic potentials with Dirac Delta 
' (— I '' couplings. The Dirac Delta coupling model has the advantage that it can be exactly solved jl3-[l3] if the uncoupled 
O |. diabatic potential has an exact solution. 

c : 

. OUR MODEL: 

^3 ■ 

CT^- . . . . 

I ! ■ We consider two diabatic curves, crossing each other. There is a coupling between the two curves, which cause 

I transitions from one curve to another. This transition would occur in the vicinity of the crossing point. In particular, 

' ^ it will occur in a narrow range of x, given by 

Vi{x)-V2ix)^Vi2{Xc), (1) 

o ■ 

' where x denotes the nuclear coordinate and Xc is the crossing point. Vi and V2 are determined by the shape of the 

10 ■ diabatic curves and V12 represents the coupling between them. Therefore it is interesting to analyse a model where 

' coupling is localized in space near Xc- Thus we put 

Vi2{x)=Kod{x~Xc), (2) 

where Kq is a constant. 



, Three simple examples 



We start with a particle moving on any of the two diabatic curves and the problem is to calculate the probability 
that the particle will still be on that diabatic curve after a time t. We write the probability amplitude for the particle 

as 

where ^i{x,t) and ip2{x,t) are the probability amplitude for the two states. The Hamiltonian is given by 

V2l{x) H22{X) ) ' 

where Hi{x), H2{x) and ^^12(2;) are defined by 

Hii{x)^^^£, + V^{x), (5) 

H22{x) = -^-l^ + V2{x) and 
Vi2{x) =V2i{x) = Ka8{x - X,). 
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The above Vi{x) and V2{x) are determined by the shape of that diabatic curve. V{x) is a couphng function which we 
assume to be a Dirac delta function. The time-independent Schrodinger equation is writen in the matrix form 



Hiiix) Vi2{x) \ f Mx) 

V2lix) H22{X) i \i>2{x) 



E 



'ii)i{x) 



This is equivalent to 



H\\{x)il:i{x) + Kq5{x — Xc)'ij^2{x) = Eijji{x) and 
Kod{x - Xc)ipi{x) + H22{x)ii2{x) = Eip2{x). 



(6) 



(7) 



Integrating the above two equations from Xc — r] to Xc + r] (where — )• 0) we get we get the following two boundary 
conditions 



2m 



dipi {x) 



dx 



+ Koip2{xc) = and 



(8) 



2m 



dip2{x) 



dx 



Xc—T) 



Also we have two more boundary condition 



'4}i{xc- rj) ^ ijiixc^ rf) and 
^2{xc -v) = ip2{xc + v)- 



(9) 



Using the above four boundary conditions we derive the transition probability from one diabatic potential to the 
other, in the case of coupling between (a) two constant potentials, (b) two linear potentials and (c) two exponential 
potential. 

Exact analytical solution for constant potential case: 

In region 1 {x < Xc), the time-independent Schrijdinger equation for the first potential is given by 

1 d^Tpiix) 



2m dx^ 

The above equation has the following solution 



(10) 



(11) 



where ki = ^ In region 2 {x > Xc), the time-independent Schrodinger equation for the first potential is given 

by 



1 d^ipiix) 



= Ei;i{x)i{x) 



2m dx"^ 

Physically acceptable solution of the above equation is given by 

V'i(a;) =Ce*'=^^. 

In region 1 {x < Xc), the time independent Schrodinger equation for the second potential is given by 

1 dH^2{x) 



2m dx'^ 



+ V2MX) = E^x). 



(12) 



(13) 



(14) 



Physically acceptable solution is given by 



ip2{x) = De 



-ik2X 



(15) 
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FIG. 1: The plot of transition probability from one constant potential to another constant potential as a function of energy of 
incident particle (7^0=1.0). 



where ^2 = \l — V2). In region 2 (a; > Xc), the time-independent Schrbdinger equation for the second potential 
is given by 



1 92^2(2;) 



+ V2i'2{x) = E'ljj2ix). 



2m dx^ 

Physically acceptable solution is given by 

Here we put = 0. Now using four boundary conditions we calculate 



and 



So, the transition probability is given by 



F 


2 




A 






D 


2 




A 




{kik2fi^ + m-^Kl) 



(16) 



(17) 



(18) 



(19) 



T = 2- 



imKQkih 



{kik2h^ + m^K§) 



(20) 



In our numerical calculation we use atomic units so that Ti — 1. In the atomic units, we set Vi{x) = 0, V2{x) = 5, 
Kq = 1.0 and m = 1.0. The result of our calculation is shown in Fig. [TJ 
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FIG. 2: Schematic diagram of the two state problem, where one linear potential is coupled to another linear potential in the 
diabatic representation. 



Exact analytical solution for Linear potential case: 



The time- independent Schrodinger equation for the case where a linear potential coupled to another linear potential 
through a Dirac delta interaction is given below (see Fig. [2]). 



■01 {x) 



= E 



ipi{x) 

^2{x) J ■ 



The Eq. (PT|) can be split into two equations 

1 



+ Pix ) ipi{x) + KoS{x)ijj2{x) = E'il)i{x) and 
P2x \ %l:2{,x) + Kf)5{x)ipi{x) =£'02 (a;). 



(21) 



(22) 



2m dx^ 

2m dx^ 

In our calculation, we use pi = P2 = 1- The Time-independent Schrodinger equation for the first diabatic potential is 
given below. 



In region 1 ( x < Xc) the physically acceptable solution is given below 

tpi{x) = {A + B)A, h^-E + x)] +i{A-B)A, h^-E + x) 



(23) 



(24) 
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Here Ai[z] and Bi[z] represent the Airy functions. In the above expression A denotes the probabihty amphtude for 
motion along the positive direction and B denotes the probabihty amphtude for motion along the negative direction 
(20j . The physically acceptable solution in region 2 ( x > Xc), is given by 



In this region the net flux is zero [20| . 

The time-independent Schrb dinger equation for the second diabatic potential is given below 



(25) 



1 

2m dx^ 

In regionl ( a; < Xc) the physically acceptable solution is 

Mx)^DA, \2^-E-x) . 
In this region the net flux is zero [20| . In region 2 ( a; > Xc) the physically acceptable solution is 



iP2{x) = F (^A, [2ii-E-x) 



iB,, 



2^^{~E~x) 



(26) 



(27) 



(28) 



Using the four boundary conditions mentioned above (here we put Xc = 0), we have derived an analytical expression 
for the transition probability from one diabatic potential to the other diabatic potential and the final expression is 
given below 



T = 16 X 2 



2/3 



where 



and 



N = Ai 



-2^'^E 



-A^ 



-2^/^E Bi -2^/^E 



+ Ai 



-2^/^E Bi' -2^/^E 



D 



-AAi 



^2^/^E 



-SiAi 

4 



'-2'/^E 


3 

Bi 


'~2"^E 


+ 2^1^ Ai! 


'~2^I^E 


2 

Bi 


'^2^I^E 



(29) 



(30) 



(31) 



22^/3^^ 



-2^I^E 

2 



Ai' 



^2^'^E 



Bi 



-2^/^E 



Bi' 



^2^/^E 



2^/^E 



-ABi 



2^/^E +2'^/^Bi' ~2^/^E 



In our numerical calculation we set pi = 1, p2 = I7 Kq = 1 and to = 1 in atomic units. The result of our calculation 
is shown in Fig. [3l 



Exact analytical solution for exponential potential case 
We start with the time-independent Schrbdinger equation for a two state system 

: E 



\ Ko6{x-Xc) -2^^ + '^o(^' 
Eq. ([32]) can be split into the following two equations 



■01 (x) 

■02(2;) 



1 92 

'2^id^ + VQe""" ] t/Jiix) + KoS{x - Xc)Mx) = Eil;i{x) and 
+ l^oe^'^" ) Mx) + KoSix - Xc)Mx) = E^x)- 



(32) 



(33) 
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FIG. 3: The plot of transition probability from one linear potential to another linear potential, as a function of energy of 
incident particle (A'o=1.0). 
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FIG. 4: Schematic diagram of the two state problem, where one exponential potential is coupled to another exponential potential 
in the diabatic representation. 
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In our calculation we took Vq = 1.0 and a = 1. Using the time-independent Schriidinger equation the first diabatic 
potential is given below 

^ +l/oe"^Vi(^) (34) 



2m dx^ 

In the region x < Xc, the solution of the above equation is given below 

V'l(2^)=^^(2.y2B)(2^^) +BI^_^^^pV2^). (35) 

Here /„ (z) represent the modified Bessel function of the first kind. In the above expression A denotes the probability 
amplitude for motion along the positive direction and B denotes the probability amplitude for motion along the 
negative direction ^20;]. In the region, where x > Xc, the physically acceptable solution is 

M^)^CK^_^^^pV2^-). (36) 

Here Kn{z) represent the modified Bessel function of the second kind. In this region the net flux is zero [l8l |. 
For the second diabatic potential, the time-independent Schrodinger equation is given below 

^ V^e-^AM^)=EMx). (37) 



2m dx^ 

In the region, where x > Xc^ the physically acceptable solution is 



V'2(x)=F/(^2^v^)(2V2^). (38) 
In the region, where x < Xc^ the physically acceptable solution is 

M^)^DKf^_^.^^pV2^). (39) 

In this region the net flux is zero [l8| . Using four boundary conditions and putting Xc = 0, as mentioned before, we 
derive an expression for transition probability from one exponential potential to the other exponential potential. In 
our numerical calculation we use atomic unit. The value of Kq = 0.1 and m = 1. The result of our calculation is 
shown in Fig. [5] 

FORMULATION OF GENERAL SOLUTION FOR TWO STATE SCATTERING PROBLEM USING OUR 

MODEL 

We start with the time-independent Schrodinger equation for a two state system 

[ V2i{x) H22{x) )[Mx)J~ \ M^) ) ■ ^ ^ 

This equation can be written in the following form 

Mx) = Vi2{x)-^[E-H^i{x)]Mx) and (41) 
= [E- H22{x)r^ V2i{x)Mx)- 

Eliminating ip2{x) from the above two equations we get 

[E - Hii{x)] Mx) - Vi2{x) [E - H22{x)r^ V2i{x)Mx) = 0. (42) 

The above equations simplify considerably if ^^12(2;) and V2i{x) are Dirac Delta function at Xc, which in operator 
notation may be written as = KqS = Kq \xc) {xc\- The above equation now become 

[[E - Hii{x)] - K^\x,){x,\ [E - H22{x)Y^ ^i{x) - 0. (43) 

This may be written as 

[i7ii(x) -I- Kl5{x ~ x,)G°(x„ X,; E)] V'i(a;) = Ei,^{x). (44) 
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FIG. 5: The plot of transition probability from one exponential potential to the other as a function of energy of incident 
particle (_K'o=0.1). 



The above equation can be written in the following form 

[E - Hn] ^Pi (x) = K^Six - x,)G°{x,, x,- E)^^ (x), 



(45) 



where the right hand side is considered as an inhomogeneous term. The general solution of this equation can be 
written as 



V'i(x) = i:a{x) + KlGl{x,Xc]E)G'i{xc,Xc;E)i;i{xc), 

where ipo(x) is a solution of the homogeneous equation 

[E Hii]iPfi{x) = Q and 
[E ~ Hi{\Gl{x,x';E) ^5{x-x'). 

In the above equation 



ipi{x) — J dx'G{x,x' ; E)ipo{x') and 

— oo 

oo 

M^)^ I dx'G°{x,x';E)Mx')- 



So Eq. can be written as 



/ dx'G{x,x';E)Mx') ^ I dx'G'i{x,x'; E)i;o{x') 

— OO — OO 

OO 

+KiG°{x,Xc;E)G''^ix,,Xc;E) J dx'Gixc,x';E)M^')- 



(46) 
(47) 



(48) 



(49) 
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The solution in terms of Green's function is as follows 

G{x, x'; E) = G\{x, x'; E) + KIgI{x, Xc\ E)GI{xc, Xc, E)G{xc,x'; E). (50) 
In the above equation we put x = x^ 

G{x„ x'; E) = G°(xe, x'; E) + K^GUxc, x^, E)G°{xe, x^] E)G{xe, x'; E). (51) 
After simplification, we get 



Gix x'-E)- G°^ixc,x';E) 

^ " ' ' 1 - K^Glix,, Xe; i?)GO(x„ x,; E) ' 



so that 



(52) 



G{x x'-E) = G?(x x'-E) + -^o^i(''^' E)G2{xc, Xc, E)Gi{xc, x'-, E) ^^^^ 

1 — KqG^{Xcj Xc', E)Gl{xc,Xc;E) ■ ' ' 

In the above expression, we have the Green's function for two state scattering problem using delta function coupling 
model. Using the expression of G(x, x'; E) one can calculate wave function and from the wave function one can easily 
calculate the transition probability from one diabatic potential to the other. 



Three channel scattering problem using our model 

We start with a particle moving on any of the three diabatic curves and the problem is to calculate the probability 
that the particle will still be in that diabatic curve after a time t. We write the probability amplitude for the particle 
as 




*(x) = M^) \ ■ (54) 



Where tpi{x,t), tp2{x,t) and ^'3(2;) are the probability amplitude for the three states. The Hamiltonian matrix of this 
system is given by 

Hnix) Vi2{x) Vi3{x) 
H= I V2i{x) H22{x) I , (55) 
V3iix) Hssix) 

where Hi{x), H2{x), H^ix), Vi2{x), V2i{x), V3i(x) and Vi3(x) are defined by 

^-^^) = -i^ + ^^^^)' 
^12(2;) = V2i{x) = KqS{x - Xc) and 

Vl3{x) = V3l{x)=KoS{x-Xc). 



In the above equation Vi(x), V2(x) and V3(x) are determined by the shape of that diabatic curve. The time- 
independent Schrodinger equation for this problem is given by 

Hn{x) KoS{x-Xc) Ko5{x-Xc)\ f Mx)\ fMx)\ 
Ko6{x-Xc) H22{x) \ \ i^2{x) \ = E \ i^2{x) \ . (57) 

Ko5{x - Xc) Hs3{x) ) \ V3(x) / V V'3(a;) / 
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This matrix representation is equivalent to the following three equations 



Hii{x)ipi{x) + KqS{x - Xc)ll^2{x) + KqS{x - Xc)i^3{x) = Eil)i{x), 
Ko5{x - Xc)iIji{x) + H22{x)ip2{x) = E'ip2ix) and 

Ko6{x - Xc)tpi{x) + H33{x)lp3{x) = El})3{x). 



(58) 



Integrating the above three equations from a;c — ?7 to aic + ?7 (where r/ ^ 0) we get the following three boundary 
conditions 



2m 



#1 {x) 



dx 



2m 



+ KqiI)2{Xc) + KolpsiXc) = 

dip2{x) 



(59) 



dx 

d^p3{x 
2m dx 



+ KqiPi{xc) = and 

Xc—rj 

- Xc+ri 

+ KoMxc) = 0. 

Xc—V 



Also we have three more boundary conditions 



■>Pi{xc-r]) =^i{xc + r]) 
4>2{xc -11)= '4>2{xc + v) and 
tpsixc-v) =i'3{Xc + V)- 



(60) 



Using the above six boundary conditions we derive analytical expressions for transition probabiUty from one diabatic 
potential to the other in the case of coupling between (a) three constant potentials, (b) three linear potentials and 
(c) three exponential potentials. 



FORMULATION OF GENERAL SOLUTION FOR MULTI-CHANNEL SCATTERING PROBLEM USING 

OUR MODEL 

We start with the time-independent Schrijdinger equation for a three state system, given by 

Hn{x) Vi2{x) Vi3ix) \ ( Mx) 

V2i{x) H22{x) I I Mx) \=E \ Mx) I , (61) 

where 



V3i{x) Hs3ix) I 




i/u(x.)^-^^ + V.(x), (62) 



This above matrix equation can be written in the following form 



[Hii{x) - E] Mx) + Vi2{x)Mx) + Vi3{x)Mx) = 0, (63) 
[H22ix) - E] xP2{x) + V2iix)Mx) = and 
[H33ix) - E] + V3i{x)Mx) = 0. 



The above equation after rearranging is given below 



[E-Hn{x)] Mx) - Vi2{x)Mx) - Vi3{x)Mx) = 0, (64) 
Mx) = [E- H22{x)y^ V2i{x)Mx) and (65) 
^3 (a:) = [E- H33ix)]-' V3i{x)Mx)- (66) 
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After eliminating both ^|J2{x) and ip3{x) from Eq. (97) we get 

{[E - Hn{x)] - VMx) [E - H22{x)]-^ V2i{x) - V^ix) [E - FgsW]"' V^iix)) Mx) = 0. (67) 

The above equations arc true for any general V12, V21, V13 and V31. The above equation simplify considerably if V12, 
V13, V31 and V21 arc Dirac Delta functions, which we write in operator notation as V12 = V21 = K2S = K2\x2){x2\ 
and V13 = V^i = K3S = K3\x3){x3\. The above equation now becomes 

{[E - Hn{x)] - K^6{x - X2)Gt{x2, X2; E) - Kl6{x - X3)Gl{x3, X3; E)) i>i{x) = 0. (68) 

This may be written as 

{[E - Hi2{x)] - KlS{x - X3)Gt{x3,X3;E)) i>i{x) = 0, (69) 

where 

Hi2{x) = Hn{x) + KlS{x - X2)G!](a;2, 2^2; E). (70) 

For Hi2{x), one can find the corresponding Green's function Gi2{x,x';E) using the method as we have used in two 
state case. 

[E - Hn]4,i{x) ^ KiSix - X2)G'^{x2,X2; E)Mx), (71) 

where the right hand side is considered as an inhomogeneous term. The general solution of this equation can be 
written as 

00 

^^{x) = ^o{x) + J dxG\{x,x';E)Kl5{x' - X2)G^{x2,X2;E)^i{x'), (72) 
—00 

where ipoix) is a solution of the homogeneous equation 

{E-Hn)Mx)=0, (73) 



where 



So 



In the above expression 



{E-Hn)G'i{x,x';E) = 5{x-x'). (74) 

Mx) = Mx) + KiGUx,X2; E)G''2{x2,X2; E)M^2). (75) 

00 

Mx)= j dx'G-i2{x,x';E)Mx') and (76) 

— 00 

00 

Mx)= J dx'G°{x,x';E)Mx')- 



So Eq. (108) can be written as 

00 

/ dx'Gi2{x,x';E)M3:')= 



00 00 

j dx'G°{x,x';E)Mx') + KiG°{x,X2;E)G°{x2,X2;E) j dx'Gi2{x2,x';E)Mx')- (77) 
—00 —00 

The solution in terms of Green's function, extracted from last equation 

Gi2{x,x';E) = G^l{x,x';E) + KiG'i{x,X2;E)Gl{x2,X2;E)Gi2{x2,x';E). (78) 
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In the above equation we put a; = a;2 to get 

Gi2{x2,x'; E) = G\{X2, x'; E) + K^Gl{x2, X2; E)G\{x2, X2; E)Gi2{x2, x'; E). (79) 
So after simplification we get 



so that 



r rr r'-F\- G\{x2,x';E) 

ui2iX2, x,i^) ^ _ ^2^(0(^2^ E)GI{X2,X2; E) ■■ 



r (r r'-F\-r^(r r' ■ F\ I KIGI{x,X2; E)Gl{x2,X2; E)G\{x2,x' ; E) 

Gx2 {x,x,E)-G,{x,x,E)+ ^_KIGI{X2,X2;E)GI{X2,X2;E) 



(80) 



(81) 



Now we will incorporate the effect of third state which is coupled to first state only, i.e. we will solve Eq. (102) in 
terms of Greens function. 

r F\-r (-r V- F\ J. ^2Gi2ix, X3;E)G°{X3,X3; E)Gi2{xz,x'- E) 

1- K^Gi2[X3,X3;E)G'^{x3,X3;E) 

So after incorporating AT-th state, where all states are coupled to the first one only we will get 

{{E-H^N_^{x))-K%6ix-XN)G%{xN,XN;E))'tPi{x)=0, (83) 

where 

N-l 

H^N-i]{x) = Hn{x)) + J2 KlS{x - Xn)Gl{xn,Xn;E). (84) 

n=2 

In this case, one can find the Green's function G\n{x,x';E) using the method as we have aheady used. 



„ , / ^ / / , KIGi[n^i\{x,xn]E)G%{xn,xn\E)Gi{n-i]{xn,x';E) 

Gin{x,x;E) = Gi[N_i]{x,x ■,E) + '—^ „o / • 

^ ' I- K^Gi[N_Vi{xN-,XN\E)G^j^{xN,XN:,E) 



(85) 



In this case wc can calculate Gin{x, xq: E) if we know Gi(Ar_i)(a:;. .tq; -E). So using Gin{x, xq: E), one can calculate 
wave function explicitely and from the wave function one can easily calculate transition probability from one diabatic 
potential to all other diabatic potentials. 



CONCLUSIONS: 



Starting from few simple examples, we have proposed a general method for finding the exact analytical solution for 
the two state quantum scattering problem in presence of a delta function coupling. Our solution is quite general 
and is valid for any potential. We have also extended our model to deal with general one dimensional multi-channel 
scattering problems. The same procedure is also applicable to the case where S" is a non-local operator and may 
be represented by 5 = |/) Kq {g\, where / and g are arbitrary acceptable functions. Choosing both of them to be 
Gaussian should be an improvement over the delta function coupling model. S may even be a linear combination of 
such operators. 
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